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We consider the role of quantum effects in the transfer of hyrogen-like species in enzyme-catalysed 
reactions. This study is stimulated by claims that the observed magnitude and temperature depen- 
dence of kinetic isotope effects imply that quantum tunneling below the energy barrier associated 
with the transition state significantly enhances the reaction rate in many enzymes. We use a path 
integral approach which provides a general framework to understand tunneling in a quantum sys- 
tem which interacts with an environment at non-zero temperature. Here the quantum system is 
the active site of the enzyme and the environment is the surrounding protein and water. Tunneling 
well below the barrier only occurs for temperatures less than a temperature To which is determined 
by the curvature of potential energy surface near the top of the barrier. We argue that for most 
enzymes this temperature is less than room temperature. For physically reasonable parameters 
quantum transition state theory gives a quantitative description of the temperature dependence 
and magnitude of kinetic isotope effects for two classes of enzymes which have been claimed to 
exhibit signatures of quantum tunneling. The only quantum effects are those associated with the 

O transition state, both reflection at the barrier top and tunneling just below the barrier. We estab- 
lish that the friction due to the environment is weak and only slightly modifies the reaction rate. 
£ — | Furthermore, at room temperature and for typical energy barriers environmental degrees of freedom 

with frequencies much less than fOOO cm -1 do not have a significant effect on quantum corrections 
to the reaction rate. 
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The possible existence and importance of quantum effects in biomolecular systems is intriguing and controversial. 
Whether quantum effects such as superposition, interference, tunneling, or entanglement are crucial to the function 
of specific biomolecules is receiving increasing attention 0, 0, H, 0, [f| H, 0, H, H E3, HU • One might expect most 
quantum effects to be destroyed by decoherence[H, [l3[ because biomolecules interact strongly with their "hot and 
wet" environment, i.e., they function at room temperature in a highly polar solvent, water. Arguably, the most 
well-established case of a quantum effect being crucial for biomolecular function is arguably electron tunneling in 
proteins [HI . Furthermore, it has been argued that by evolution electron transfer proteins vary and are selected based 
on tunneling parameters [l4j . The role of tunneling in other biomolecular systems has also been examined [15]. For 
example, in myoglobin it has been found that the reaction rate for binding of carbon monoxide becomes independent 
of temperature below 80 K, due to the presence of quantum tunneling (la. [l7|. 

Over the past two decades the possibility of quantum tunneling of protons in enzymes has attracted considerable 
attention [TJ [l8|, [2(| HH . The large kinetic isotope effects and their temperature dependence are inconsistent 
with semi-classical transition state theory, whereby the chemical reaction occurs via thermal activation over an energy 
barrier. These discrepancies have been interpreted as evidence for the presence of tunneling (Til. I22I l23l. l24l l25l. [2rj |27|. 
(H, US US HH HI, [H[ [13, HH . However, it should be stressed that this evidence is rather indirect, being based on the 
values of fitting parameters for Arrhenius plots for the temperature dependence of the reaction rate, where the absolute 
temperature only varies by about ten per cent. In contrast, for chemical reactions involving much simpler organic 
molecules, such as benzoic acid[36| or hvdroxvmethvlene [37) much more definitive signatures of proton tunneling have 
been observed. These include a temperature independent rate at low temperatures and tunnel splitting of the ground 
state energy [H, [39j] . 

Key questions that need to be answered include: 

Can some of the experimental results be explained without invoking tunneling? 

To what extent is it necessary to go beyond semi-classical transition state theory to explain the observed kinetic 
isotope effects? 

If tunneling does occur, is it important for the function of the enzyme? 
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Have enzymes evolved in a manner that enhances the contribution of tunneling? 

There are currently a wide range of views on the answers to these questions. For example, a review in Science 
states that, "the entire and sole source of the catalytic power of enzymes is due to the lowering of the free energy 
of activation and any increase in the generalized transmission co-efficient, as compared to that of the uncatalyzed 
reaction" [40(. Villa and Warshel state that, "the most important contribution to catalysis comes from the reduction 
of the activation free energy by electrostatic effects ... the popular proposal that enzymes catalyze reactions by 
special dynamical effects is not supported by a consistent simulation study ... the interpretation of recent experiments 
as evidence for dynamical contributions to catalysis is unjustified." [4l| . In contrast, Klinman et al. state that, 
" Our present findings on hydrogen transfer under physiological conditions cannot be explained without invoking both 
quantum mechanics and enzyme dynamics." [30| . Furthermore, Klinman and Kohen proposed that, "The optimization 
of enzyme catalysis may entail the evolutionary implementation of chemical strategies that increase the probability of 
tunneling and thereby accelerate the reaction rate." [12] However, Doll, Bender, and Finke[43|,|44[ synthesized artificial 
catalysts which performed the same chemistry in solution (i.e., in the absence of the enzyme) and exhibited the same 
kinetic isotopic effects. In a paper that focused on simulations Schwartz et al. [45] express a similar view to Klinman's, 
"The action of the enzyme in speeding the chemical reaction, however, is postulated to be intimately connected to the 
directed vibrational motion identified in this paper. Thus, it appears that evolution has designed the protein matrix 
of an enzyme not just to hold substrates or stabilize transition state formation, but rather to channel energy in a 
specific chemically relevant direction." . 

Over the past five years several reviews of different theoretical approaches to this problem have appeared [Til HH, |46|. 
H3, H^|. Most theoretical work makes two particular assumptions which may be debatable, (i) that the proton transfer 
process is adiabatic, and (ii) that a sing le reaction co-ordinate is adequate. For a detailed discussion of these issues we 
refer to a recent review by Marcus, [lg We also note that for non-enzymatic reactions, the first assumption has been 
brought into question and an alternative non-adiabatic picture (analogous to electron transfer) has been proposed. [4!| 
The non-adiabatic proton transfer theory has been applied to enzymes. [13, HH, [H, [H[ Then the only way for the 
proton to move from the reactant to product state is via tunneling. The activation energy is then associated with the 
reorganisation of the environment rather than that of the transition state. Siebrand and Smedarchina[54j considered 
such a approach to explain how some enzymes have a large KIE that is weakly temperature dependent. 

One approach to examine the role of quantum effects in complex biological molecules is to use quantum mechanical 
molecular mechanics (QM-MM) simulations. In this approach atoms which are directly involved in the reaction are 
treated quantum mechanically while the rest of the enzyme is treated classically. This approach has been applied to 
a number of different enzymes |63l [H, [65| . Schwartz and coworkers have as their starting 



point a Hamiltonian similar to the one used here [45|, 47l l66l|. They used classical molecular dynamic techniques to 
simulate specific reactions and extract the spectral density. 

Using a low-energy effective Hamiltonian model such as the Caldeira-Leggett Hamiltonian [67j to capture the essen- 
tial physics of the relevant process offers a complimentary approach to QM-MM simulations. It has the advantage that 
quantum effects and the role of the environment can be treated more rigorously, via path integral methods (l2l. [8^. 
Furthermore, the dependence of behavior on the key physical parameters such as the curvature of the potential energy 
surface near the transition state can be elucidated. 

In this paper we establish the following points using a path integral approach, (i) Tunneling well below the barrier 
only occurs for temperatures less than a temperature Tq which is largely determined by the curvature of the top of 
the barrier. We argue that for most enzymes this temperature is less than room temperature, (ii) For physically 
reasonable parameters quantum transition state theory gives a quantitative description of the temperature dependence 
and magnitude of kinetic isotope effects for two classes of enzymes which have been claimed to exhibit signatures of 
quantum tunneling. The only quantum effects are those associated with the transition state, both reflection at the 
barrier top and tunneling just below the barrier, (iii) The friction on the proton due to the environment is weak and 
only slightly modifies the reaction rate, (iv) At room temperature environmental degrees of freedom with frequencies 
much less than 1000 cm -1 do not have a significant effect on quantum corrections to the reaction rate. 

II. BACKGROUND 
A. Kinetic isotope effects 

The rate coefficient fc^, for a chemical reaction involving a species L at temperature T can be written in the 
Arrhcnius form 



k L = A L exp{-E L /(k B T)), 



(1) 
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where El denotes the activation energy for the reaction and Al is the prefactor. The two quantities Al and El are 
generally referred to as the Arrhenius parameters. 

The reactions that we will be interested in all involve breaking or forming bonds which contain hydrogen species 
(protons, deuterium, tritium, hydrogen atoms, and hydrogen anions). We will only be considering the primary kinetic 
isotope effects for systems where the hydrogen transfer step is rate limiting. The kinetic isotope effect (KIE) is 
generally expressed as the ratio of rate constants kn/kx or kn /&D, where the superscripts P, D and T denote the 
reactions in which a proton, deuterium, and tritium are being transferred, respectively. 



B. Semi-classical Transition State Theory 



Consider the one-dimensional potential energy shown in Figure 1. A reaction from A to C proceeds via a transition 
state at B. The following expression for the rate coefficient is widely used [68, 69], 

— ^— J 7^- exp (~E b /k B T) = k k T sT, (2) 

where Za is the quantum partition function of the metastable state A, and Z$ is the quantum partition function of 
the activated complex B. The parameter k was originally introduced to account for the fact that some trajectories 
may recross the transition state and return to the reactant state and is also set used to include the effects of quantum 
tunneling. An ad-hoc semi-classical transition state theory the activation energy is replaced by Ef, — Jtwq/2 which 
corrects for the effect of the quantum-zero point motion in the reactant well. This leads to kinetic isotope effects 
because luq depends on the mass of the hydrogenic species being transferred. 

Kim and Kreevov[70j gave three criteria that are widely considered to be signatures of quantum tunneling in 
hydrogen-transfer reactions: (i) a deuterium kinetic isotope effect kn /ko significantly larger than 6.4 at 20 C (or 8.9 
if secondary isotope effects are included); (ii) an activation energy difference, Eh — Ed larger than 5.0 kJ/mol; and 
(iii) a ratio of prefactors, Ah/Ad less than 0.7. However, it should be noted that Kim and Kreevoy also stated, "it 
appears that completely unambiguous experimental proof that tunneling occurs at ~ 300 K would be impossible to 
obtain, although tunneling clearly becomes evident at much lower temperatures." 

Following Bell [7lj], Kohen et glff^ state that if tunneling is not significant prefactor ratios should fall within the 
range of 0.3 < A H /A T < 1.7 and 0.5 < A D /A T < 1.4. 



C. Heuristic justification for Quantum Transition State Theory 

Here we reproduce simple arguments described by Weiss (l2|. A particle in thermodynamic equilibrium in the 
reactant well A is in a metastable state and so we can think of each quantum state that contributes to the system 
partition function Z has an imaginary part, i.e. e„ = E n + ihT n , with E n 3> hT n , and 

Z = £exp(- e „/fe s T) = Z x +iZ 2 ~ £ exp(-E n /k B T) - i £ ^ exp(-£„/fe B T). (3) 

n n n 

The total decay rate out of the reactant well is then 

fc = ^£r„ex P (-^ B T) = ^§. (4) 

If the motion in the reactant well A is described by a single harmonic oscillator in thermal equilibrium, with 
frequency u)q, the partition function is, 

(5) 



sinh 



hue* 



If the barrier is an inverted parabola then in this partition function we can replace loq with iu>b to obtain 

Z* = -A ^ exp (-E b /k B T). (6) 



sin 
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FIG. 1: Potential energy as a function of the reaction coordinate, x, with the metastable reaction state at A, the transition 
state at B, and the final product state at C. Escape occurs via the forward rate k and Eb is the corresponding activation 
energy. The angular frequency of oscillations about the reactant state is loo, which depends on the curvature of the potential 
energy surface at the local minimum (x = x a ) and the mass of the particle. Similarly the barrier frequency u>b, depends on 
the curvature of the potential energy at the local maximum (x = Xb) and the mass of the particle [68| ] . This paper addresses 
the question as to what extent the reaction A — > C can proceed via quantum tunneling below the barrier for enzyme catalysed 
hydrogen transfer reactions. 



Assuming no quantum coherence between the bottom and the top of the barrier the total partition function is then 
Z = Za + iZ*. Substituting this in (0} then gives 

, , sinh (^p) 

k = ^ ^1 exp (-E b /k B T) (7) 
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An important limitation of this expression is that it is only well defined for temperatures T > Tq where 

To = (8) 

The expression (J7|) was actually derived by Wigner in 1932[73j , using an expression for the energy dependence of the 
transmission probability through a parabolic barrier. This takes into account the fact that in quantum mechanics 
a particle with energy E > E b has less than unit probability of transmission (i.e. above barrier reflection occurs). 
Wigner also assumed that T ^> Tq. Bell[7_l( also derives and discusses this expression. Note that in the limit that 
hu>b <C 2ksT <C Hloq, ((7]) reduces to the semi-classical expression 

t rrt 

k = ^-exp-(E b -hu> /2)/k B T). (9) 
h 

In Section IlII El we use path integral methods to give a rigorous derivation of Wigner's expression which will also 
elucidate its range of validity and the physical significance of the temperature scale To. In Section IVlII we show that 
this expression can be used to give a quantitative description of the magnitude and temperature dependence of kinetic 
isotope effects in several important classes of enzymes. In Section [V] we will see the fact that the rate diverges as T 
approaches T is an artefact of treating the potential barrier as parabolic. 



D. Arrenhius parameters for enzymes are inconsistent with Semi-classical Transition State Theory 



Table U summarises the experimentally determined kinetic parameters for a number of enzymes. The evidence 
for tunneling generally comes from examining the prefactor ratios and also the difference in activation energy for 
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different isotopes. Specifically, when these quantities lie outside the bounds proposed by Bell[7l| it is usually claimed 
that tunneling occurs. As a reference the activation energy for the different reactions have also been included. The 
Table also includes some data for other organic reactions which have kinetic parameters that are inconsistent with 
semi-classical transition state theory. 



TABLE I: Deviations of the Arrhenius parameters for the kinetic isotopic effects of hydrogen transfer reactions in a range of 
enzymes from the predictions of semi-classical Transition State Theory. A number of hydrogen transfer reactions involving 
small organic molecules also exhibit parameters that fall outside the semi-classical limits. 
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III. RATE THEORY 

A. The Caldeira-Leggett Model Hamiltonian describes a quantum system interacting with its environment 



Consider a system which consists of a single particle of mass M described by one degree of freedom and coupled 
to a large environment which can be represented by a bath of harmonic oscillators. This is equivalent to representing 
some arbitrary environment in terms of its normal modes. The interaction of the degree of freedom with each of the 
bath modes is inversely proportional to the volume of the bath. Hence for a macroscopic environment the coupling 
to each of the individual modes will be weak [T3|. The Hamiltonian can be represented as 



H = 1~Ls + 'He + Ttz, 

where 

Hs = P 2 /2M + V(x) 

is the Hamiltonian associated with the reaction coordinate x. The Hamiltonian of the environment is given by 

N ' ,2 



+ m a ulql , . 



(10) 
(11) 

(12) 



and describes N harmonic oscillators where where m a and u a are the mass and frequency of the ath oscillator. The 
interaction Hamiltonian is 



Hx = -} j F a (x)q a . 



(13) 



In our case we will require that the interaction is separable, i.e., F a {x) = C a F(x), and that the dissipation is strictly 
linear, i.e., F(x) = x. This describes state independent dissipation, i.e., the magnitude of the friction is the same at all 
points along the reaction coordinate. In that case the complete Hamiltonian is the Caldeira-Leggett Hamiltonian [67| : 



n 



r> 2 K 1 N 



pi 



m a L0 o 



where we have introduced an effective curvature induced by the environment and defined by 

N ^2 



We also define the effective potential 



K, 



E 

a=l 



m a oj : 



2 ' 



U{x) = V{x) + ^fx 2 . 



(14) 



(15) 



(16) 



B. The spectral density describes the frequency-dependent friction due to the environment 



The role that the environment plays can be embodied in a single function. It depends on how the coupling strength 
to each oscillator mode changes with the frequency of the oscillator. This can be expressed in terms of the memory 
friction kernel, ^(t) defined aspjj, 



i N r 2 

\/l < J m _ /,j z 



a—l " 



(17) 



The Laplace transform of the memory friction kernel, 7(2), is what determines the effect of the environment on the 
reaction rate, 



1 N r 2 



M ^— ' m n u). 

a — l 



(18) 
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The spectral density, J(u>) is defined as 

a 

and is an alternative means of characterising the coupling to the environment. The spectral density and the Laplace 
transform of the memory friction kernel are related by the identity. [12l] 

* 2 f°° , 1 

Mtt J uj lo 2 + z 2 K ' 

From this we can obtain an upper bound for the friction kernel 

2 f°° , JUS) 2K e , s 

l(z) < —rr- / du——!- = — - 21 

where K e is the curvature (| 1 5[) . This expression will be used in Section I VII to estimate the magnitude of the friction. 

The simplest kind of dissipation is memoryless friction, j(z) = 7 or J(u>) — Mjuj. However, there is always some 
microscopic memory which sets the time scale on which inertial effects in the bath are significant. The simplest form 
of damping kernel that captures this is the Drude regularisation. [l2]| 



7(z) = 



1 + z/uj d 



(22) 



1 + LU 2 /UJ 2 D ' 

We can model a peak in the spectral density at a frequency tu r by 



Re 7 M = # = , 2 7r( " r)2 . ^ (23) 



which has a value of j r at the peak which has a width T. The corresponding friction kernel is 

= 24 , 7r f^ r - (24) 

C. Classical Kramers Theory defines an effective barrier frequency 

We first review results for the classical limit of the Hamiltonian. After averaging over all the environmental variables 
one finds from Eq. (fT4|) the generalised Langevin equation 

Mx+^-+M f ~/(t-s)x(s)ds = £(t), (25) 

where £(t) is the random force the particle experiences and 7(4) is the friction kernel that describes the dissipative 
interaction with the environment. When the total system is prepared initially in thermal equilibrium, the random 
force becomes a stationary Gaussian noise of vanishing mean, i.e. = 0. The classical fluctuation-dissipation 

theorem gives 

mm) - Mk B T>y(t), (26) 

where the friction kernel j(t) is defined in Eq. (|17[) . From this one can then perform a normal mode analysis to 
evaluate the partition functions entering the transition rate expression [l2j • Close to the bottom of the potential well 

U(x)n±MLjl(x-x ) 2 , (27) 

and at the barrier top 

U{x) &E b - \Mu) 2 b (x - x b ) 2 . (28) 



8 



This yields the classical rate which incorporates a dissipative interaction with the environment [68j 

fcd = — ^exp(-£ b /fc s T), (29) 
and here the effective barrier frequency \i is the solution of the equation, 



, = V f ^+^- 1 T- (30) 

In this framework the activation energy has no mass dependence and hence remains unchanged by an isotopic 
substitution. The only quantities that are altered are too, 0Jb and /i, which all appear in the prefactor. The particle 
mass term appearing in luq and ujb cancel each other out so that the only mass dependence lies in the effective 
barrier frequency. The entire KIE comes from the effective barrier frequency. Eq. (|30p then gives the bounds on the 
KIE. In the case where one is comparing the rate of a reaction where protium is transferred with a reaction where 
tritium is transferred this predicts that 1 < kn/kx < 1.7. Similarly for protium and deuterium 1 < ku/kjj < 1.4. 
Experimentally the KIE does often depend on the temperature and in the systems of interest often falls outside of the 
former bounds (See Table HI . The inconsistency between the KIEs predicted by classical Kramers theory and those 
measured by experiment shows that this classical description does not capture the relevant physics. 

D. Path integral representation of Quantum Kramers Theory 

The general problem of quantum tunneling at non-zero temperature in the presence of an environment can be 
treated using complex-time path integrals [l2l l68l [85| . Consider the partition function Z = Tr{exp(— H/ksT)}, 
where Ti. denotes the full Hamiltonian operator corresponding to the system plus environment. This quantity can be 
expressed in the form of a functional path integral over the tunneling coordinate x(t) [6§j |. 

Z = J Vx{T)exp{~S E [x(T)}/h}, (31) 

where r = it is a real variable. This integral sums over all paths x(t) that satisfy the periodic boundary condition 

x(t = -6/2) = x(t = 9/2). (32) 
with period 9 determined by the temperature, 

9 = h/{k B T). (33) 

After integrating over the bath modes the effective Euclidean action takes the form, 

r e/2 M 1 r e/2 r e/2 

S E [x]= d T{—X 2 (T) + U[x(T)}} + - dr dr'C(T-T') X (T)x(T'). (34) 

J-8/2 * * J-8/2 J -8/2 

The influence kernel £(r) is periodic in imaginary time with period 9. It is related to i(z), the Laplace transform of 
the memory friction (see Eq. (JT5J)) and can be represented in terms of a Fourier series as [86| 

M °° 

C(t) = y Kl7(KI)exp(w„T) (35) 

n— — 00 

satisfying 

' C(r)dr = 0, (36) 

-8/2 

and where v n are the Matsubara frequencies for bosons, 

v n = n2nkBT/h. (37) 
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For a metastable potential the partition function has an imaginary part which can be related to the escape rate 
from the potential [12J. The dominant contributions to the partition function, and indirectly the rate expression, come 
from the vicinity of paths in which the action (|34|) is stationary. These paths, x e (r), satisfy the equation of motion 

Mx e (r) - dU J Xe }\ )] - f' 2 dr'Cir - r')x e (r') = 0, (38) 

OX e {T) J_ e/2 

and the periodic boundary condition (|32p . In the absence of dissipation, i.e. £(r) = 0, the evolution of x e (r) in 
imaginary time corresponds to real-time motion in the metastable inverted potential — C/(x).[13, [89| 

Because of Eq. (|36|) . Eq. (|38|) has two trivial but physically important solutions. The first where the particle 
remains at the bottom of the reactant well (x e (r) = x a ) and the other where it sits on top of the barrier (x e (r) — xj,)- 
The latter corresponds to thermal activation over the barrier top. 



E. The bounce solution describes quantum tunneling which only occurs below a temperature, To 

A non-trivial periodic solution to Eqn. (138|) (which has been dubbed the bounce or instanton solution) describes 
quantum tunneling below the barrier. [87l l88l| This solution only exists below a temperature Tq [68| . In the absence 
of dissipation, an analytic solution has been found for an inverted parabola, an Eckart potential. [89]] and a cubic 
potential. [90l] In the presence of Ohmic dissipation, an analytic solution for a cubic potential has been found for 
specific values of the dissipation. j9lj 

For temperatures T > To the period of the 0-periodic orbit is not of sufficient length to admit an oscillation of the 
particle in the classically forbidden regime. For the case of zero temperature the bounce solution gives a rate which 
is directly proportional to the tunneling probability calculated from the WKB approximation, 

T{E) = exp (-2S(E) /h) (39) 

where S(E) is the value of classical action along the imaginary time trajectory, 

r-x-2 

S(E) = V2M dx[U(x)-E} 1/2 , (40) 



where x\ and X2 are the classical turning points for energy E. 

For temperatures T > To the bounce solution does not exist and the only contribution to the path integral comes 
from the constant solution (x e (r) — Xb) where the particle sits at the barrier top. Fig. 40 in Ref. [H depicts the 
different kinetic regimes which occur for different temperatures. 

We now focus on the case of a parabolic barrier. In the absence of any dissipation (7 = 0), the crossover temperature 
has the value given by ©. In the presence of dissipation, the crossover temperature T is given by 

To ^^- = 0.23 K-^, (41) 
Zttkb cm 1 

where \x is the effective barrier frequency defined in Eq. (|30|) . This means that tunneling can only occur at room 
temperature if /i > 1300 cm -1 . In an appendix it is shown how for a Lorentzian spectral density fj, is reduced by 
friction. 



IV. QUANTUM CORRECTION FACTOR: T > T 



In the high-temperature regime for a parabolic barrier one can obtain an analytic expression for the rate constant [92j, 



fc(T) — kclCqr, 



— ( 


/ w \ 







luq + n 2 v 2 + nv^j(nv) 
-J^ —uj 2 + n 2 v 2 + nvj(nv) 



n 



The first term in the square brackets denotes the classical Kramers rate for memory friction (Eq. 
smallest Matsubara frequency, 



(42) 

Here v is the 



v = 2nkBT/h 



(43) 
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and we must have v > // where fi is the effective barrier frequency given by (|30|). We note that (|42p is proportional 
to Z$ /Za where the corresponding partition functions are for damped harmonic oscillators. [I2j 

The quantum correction to the rate expression, c qm , is encapsulated by the term inside the curly brackets. For 
T 3> To this correction factor approaches unity. Moreover, it always exceeds unity which implies that quantum effects 
always enhance the classical rate. 



A. The weak friction limit reduces to Wigner's expression 

This limit has to be treated with some care as when the friction is exceptionally weak thermal equilibrium no longer 
prevails in the reactant well. It has been shown that so long as the following condition is satisfied the assumption can 
be made that the reactant system is always in thermal equilibrium 68], 

m > w (44 ) 

However, we note that for some of the enzymes shown in Table I the activation energy is sufficiently small this 
assumption may not be justified. 

In the limit where 7^0 Eq. (|4"2"|) can be simplified such that the correction factor can be written as 



This expression will be a reasonable approximation provided that for all n = 1, 2, . 



(45) 



^ < 1. 46 

nv 

Since v ~ 1300 cm -1 at room temperature this means that any friction associated with environmental modes of much 
lower frequency may have little effect on the quantum correction factor. 

Note that the expression (|4"5|) diverges as T — > T + . This divergence turns out to be an artefact from treating the 
potential near the transition state as a perfect inverted parabola. Below we will show how a more a rigorous treatment 
shows that in this temperature regime the correction factor is always finite. Furthermore, for realistic parameters 
is a good approximation down to 1.1Tb. (See Fig. [5]). 



B. Apparent Arrhenius parameters in the weak friction limit 



Most kinetic experiments are performed over a narrow temperature range. The temperature dependence appears 
to be activated (i.e, a plot of ln(fcf,/fcr) vs. 1/T is linear over the measured temperature range) and so it is natural 
to determine the Arrhenius parameters. If the full QTST rate expression is expanded about room temperature one 
can obtain an expression for the KIE that has a simple activated temperature dependence. The subscripts L and 
denote possible combinations of the three different isotopes of hydrogen. This is typically what is done in experiments 
where the heavier isotopes are used as a reference. Combining the results from Eq. (|42]) and Eq. (|4"5]) gives an 
expression for the rate constant. The KIE is given by 




where mi is the unitless mass number of the L isotope and uj$ and are the ground state oscillation frequency and 
barrier frequency, respectively for hydrogen. For a typical C-H stretch frequency, Hujq 3> 2fcsT, at room temperature 
so the hyperbolic sine terms can be approximated as exponential functions. Over the biologically relevant temperature 
range T only varies by less than 10% and so one can expand the other temperature dependent parts of the expression 
up to linear terms in 1/T. This gives an expression for the KIE which has a simple activated temperature dependence 
with the following apparent Arrhenius parameters. 
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x exp 



m Lt sm(f3 R hu} b /y/4m L t) 
m L sm(/3 R hu; b /y/<im L ) 
(3 R hwb ( cot({3 R fiw b /y/4m L t] 



cot(/3 R hw b /y/imZ) 



E 



Lt 



E, 



+ 



2 V 

hujQ / l l 

2 \*/mZ y/ml 
fuv b ( cot(f3 R tujj b /yJAm L ' i 



(48) 



cot{l3 R frw b /^JAm L ) 



where Tr = l/hupR is the temperature around which the expansion is performed. From these expressions it is also 
possible to place bounds on the apparent Arrhcnius Parameters. When T R > Tq and mj, < m^t it can be shown that 
the Arrhenius parameters are monotonic functions of (3 R fuv b . 

On their own the expressions in Eq. (148|) may not seem to shed much light. However, evaluating these expressions 
for typical parameter values shows that typical values for ujq and uj b give kinetic parameter values that are inconsistent 
with the predictions of standard semi-classical rate theory. Moreover, Table HT1 shows that the parameter trends are 
consistent with experimentally determined Arrhenius parameters for a number of systems. In a number of systems 
where tunneling has been invoked the difference in apparent activation for the different isotopes is greater than would 
be predicted by semi-classical theories. The effective activation energies derived from QTST are quantitatively similar 
to a number of the experimental values and significantly exceed the semi-classical values. The prefactor values that 
have been observed experimentally are smaller than would be expected from semi-classical theories. 



TABLE II: The effective Arrhenius parameters calculated for the expanded QTST expression are compared with experimentally 
determined values. Both the calculated and experimentally determined parameters are inconsistent with a standard semi- 
classical analysis. This shows that there are a subset of enzymatic systems where the anomalous kinetics can be explained by 
Quantum Transition State Theory. 



Parameter Values/Experimental System 


E T - E H (kJ/mol) 


A H 1 A t 


E T - E D (kJ/mol) 


A D 1 A t 


Semiclassical limits 

luo = 3000 cm" 1 , u b = 1000 cm" 1 , Tr = 288K 
Bovine Serum Amine Oxidase [23] 
Flavoenzyme Monoamine Oxidase B [25] 


< 10.1 
16 

14.2 ± 0.7 
13 


0.3-1.7 
0.08 
0.12 ± 0.04 
0.13 


< 3.1 
3.6 

4.51 ± 0.48 
4.1 


0.5-1.4 
0.7 
0.51 ± 0.1 
0.51 



V. NON-PARABOLICITY OF THE BARRIER ONLY MATTERS AT TEMPERATURES CLOSE TO T 

In the regime where T ~ To the action associated with the bounce solution and the trivial solution at the top 
of the barrier become comparable (See Section IIII Ep . In the limit where T — * Tq the expression for the quantum 
correction given in Eq. (|42|) diverges. This divergence is a consequence of the fact that the saddle point approximation 
of the imaginary time functional integral employed in evaluating the expressions presented so far breaks down in the 
vicinity of T). In order to obtain an expression which is valid around the crossover temperature the analysis needs to 
be extended to include the effects of a non-parabolic barrier [86|. The potential is expanded around the barrier top to 
give 

V(x) =E b - \Mul{ x Xb f + £ MCfc(a p b)fc . (49) 

k=3 

In the weak friction limit (see Section IIV Ap the rate expression can be simplified to give a quantum correction 
factor which takes the form 

u) b smh(huJo /2k bT) 

CqMT = ~ sm(hw b /2k B T) (50) 
x V7r(-e)(l - e/2)reerfc(-e(l - e/2)n) exp(e 2 (l - e/2)V). 

where 



e = (T - T)/T . (51) 
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This extension is only required when the temperature is close to the crossover temperature. The extent of this 
crossover region depends on both the scaled temperature e, and the parameter k, defined as 



K = U? b 



8M 
Bk B T 



(52) 



where 



B = 4c%/3u% + 3c 4 . 



(53) 



If the barrier is high and the potential reasonably smooth then Muj\ /B is an energy which is of the order of the barrier 

However, 



height so that k is of the order of yj (Eb/ftub) ^ l[86j|- However, caution is in order because for Eb ~ 50 kJ/mol 
~ 4000 cm' 1 this condition is only weakly satisfied. Also, quantum chemistry calculations discussed in Section IVl Bl 
show highly non-parabolic barriers. 
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FIG. 2: Comparison of the approximate quantum correction factor with the exact result. The approximate expression does not 
take into account the non-parabolic character of the barrier top. For temperatures above about 1.1Tb the two expressions (|45p 
and (|50[) are in good agreement. However as To is approached the high temperature expression diverges. The parameter values 
used in generating these plots are: loq — 3000 cm -1 , = 1000 cm -1 , k(T ) = 10. This shows that we are justified in using 
our expression (|62|) to describe the kinetic isotopic effect provided that the crossover temperature To is less than about 250 K. 



VI. ESTIMATES OF MODEL PARAMETER VALUES FOR ENZYMES 
A. Frequency of Oscillations in the Reactant Well (i^o) 



The parameter ujq is defined as y/U"(xo)/M, where M is the effective mass of the hydrogen isotope being transferred 
and U"(xq) is the curvature of the potential at x = xq. Using IR spectroscopy the oscillation frequencies of the different 
chemical bond stretches can be measured experimentally [94] . Results from molecular dynamics simulation show that 
the value of ujo is similar to the value obtained from IR spectroscopy. In [58| it was established that the oscillation 
frequency of the C-H bond that participates in the reaction, ujo is 2900 cm -1 . A typical IR spectrum of an organic 
molecule would contain peaks corresponding to C-H stretch frequencies in the range 2700-3300 cm -1 [941 ]. 

The value of ujq has the potential to be influenced by interactions with atoms in the active site. This interaction 
has been known to decrease the stiffness of the bond j95|. As a result it is possible that in some of the systems we 
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are studying the binding of the substrate could bring about a change in the ground state oscillation frequency, but 
based on what happens with hydrogen bonds it is unlikely this could ever be a reduction larger than a factor of two. 

B. Barrier Frequency (u>(,) 

The barrier frequency, u>b, strongly influences to what extent quantum effects affect the kinetics of the hydrogen 
transfer process. It depends on both the curvature of the potential at the barrier and also on the mass of the hydrogen 
species being transferred. It is a parameter which would be difficult to obtain experimentally. Generally one needs 
to resort to quantum chemistry calculations to obtain an estimate of u>b- For our purposes u>b is a crucial parameter. 
Firstly, it sets the upper limit of the crossover temperature, To, around which the kinetics becomes more classical 
in nature. Secondly, in the intermediate and high temperature regime it strongly influences exactly how much the 
rate gets modified due to quantum corrections. Table IIIII shows the values of the barrier frequency that have been 
calculated using quantum chemistry for a number of different enzymes. 

TABLE III: Barrier frequency values for several enzymes estimated from quantum chemistry calculations of potential energy 
surfaces for the hydrogen transfer reaction. The maximum crossover temperature To is related to cub by Eqn. JSJ. Friction will 
tend to reduce these values. The corresponding temperatures for deuterium and tritium transfer will be about 30 and 50 per 
cent lower, respectively. SCC-DFTB denotes Self-Consistent Charge-Density Functional Tight Binding. AMI and PM3 are 
semi-empirical methods, and SRP denotes a Specific Reaction Parametrisation, against a DFT calculation using the B3LYP 
functional and the 6-31G* basis set. 



Enzyme 


Level of theory 


Wf,(cm i ) 


Max To (K) 


Ref. 


Triosephosphate Isomerase (TIM) model (gas phase) 


AM1-SRP 


1365 


315 


f55] 


TIM model (in water) 


AM1-SRP 


591 


140 


[55] 


TIM model (in protein model) 


AM1-SRP 


798 


190 


f55] 


Liver Alcohol Dehydrogenase(LADH) 


SCC-DFTB 


783 


180 


f57] 


LADH 


AMI 


1046 


240 


[58] 


LADH 


AMI 


1229 


240 


[63] 


Monoamine Oxidase 


B3LYP/6-31G* 


1054 


240 


[96] 


Monoamine Oxidase 


PM3 


1782 


410 


[96] 


Methylamine dehydrogenase (MADH) 


PM3 


2000 


460 


[60] 


MADH 


AM1-SRP 


2218 


510 


[63] 


MADH 


PM3 


2000 


460 


[97] 


Soybean Lipoxygenase 


PM3/d 


2913 


670 


[63] 


Aromatic amine dehydrogenase (tryptamine) 


PM3-SRP 


2057 


450 


[98] 



From these calculated values of ujb it is clear that in a number of the cases considered the barrier frequency is of 
the order of 1000 cm , which corresponds to a maximum crossover temperature of around -40 °C. Of the enzymes 
that appear in the table both Methylamine dehydrogenase and Soybean Lipoxygenase have a large KIE, associated 
with a prefactor ratio which is much larger than one (See Table HJ. 

It needs to be kept in mind that the values for the the barrier curvature need to be taken with some caution. Most 
of the values listed in Table [TTT1 were determined using hybrid QM/MM calculations. These techniques encounters 
some methodolog ical problems which have the potential to strongly influence the calculated barrier frequencies [99j . 

A recent study [100| performed a QM-MM study of the hydrogen transfer reaction in methylmalonlyl-CoA mutase, 
treating 44 atoms near the active site at the AMI level. They emphasize the role of corner cutting but do not give a 
value for the curvature of the barrier. The "representative tunneling energy" (the energy at which the product of the 
transmission coefficient and the Boltzmann factor is a maximum) is about 1400 cm -1 below the top of the barrier. 

Obtaining a reliable value for the barrier curvature from computational chemistry represents a major challenge. We 
note how the Table shows that the values obtained depend on the level of theory. Other factors to consider are the role 
of anharmonicity, dependence of the results on the active site geometry, and the fact that in a dynamic environment 
(protein plus water) there are actually many reaction paths. 

C. Frequency-dependent Friction due to the Solvent-Protein Environment 

The analysis in the previous Sections shows that the magnitude and frequency dependence of friction j(u>), and its 
Laplace transform, the memory friction kernel 7(2), can lead to not just quantitative but also qualitative differences 
in the reaction rate. For example, if j(uJb) > w fci then the friction can substantially reduce the temperature T (see 
Figure^ below which the "bounce" solution exists, i.e., strong enough friction can prevent tunneling from the bottom 
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of the well oceuring. The friction also causes vibrational energy relaxation and dephasing of the vibrations in the 
reactant well, jioij The former is proportional to Re j(ujq). Henc e, m easurements of the relaxation and dephasing 
rates provides a means to determine the magnitude of the friction. |l02l | 

There are believed to be several main sources of friction associated with bond deformation and breaking in 
biomolecules. The first source is the interaction of the dipole moment associated with displacement of the pro- 
ton with the fluctuating electric field associated with fluctuations in its environment . 1 1031 ] the surrounding protein 
and solvent. Recent ultrafast infrared spectroscopy experiments of the OH stretch of HOD in liqui d D 2 O have s hown 
that in bulk water this fluctuating electric field is the dominant source of vibrational dephasing. [104L Il05l . llOq The 
second source of friction is the interaction of oscillations of the reaction co-ordinate with a modulating low frequency 
mode, w hich in t urn is strongly damped by the environment. This is particularly important in hydrogen bonded 
systems. [l02l Il07| For the case, A-H- • • B, the A-H stretch is modula ted b y the A- • • B oscillations associated with the 
hydrogen bond. A third possible source of friction is anharmonicities [l08j | and Fermi resonances with other vibrational 
modes in the biomolecule. 



1. Comparison with the spectral density for biological chromophores 



The dynamics of optically active molecules (chromophores) within proteins have been studied extensively, both 
theoretically and experimentally. [6[ An optical transition from the ground to an excited electronic state is usually 
associated with a change in electric dipole moment, A/x of the order of a few Debye, which then couples to the 
electric field (reaction field R(t)) associated with the dielectric relaxation of the chromophores environment. The 
relaxation can usually be assigned to three components of the environment: the surrounding protein, water bound to 
the surface of the protein, and the bulk water surrounding the protein. The corresponding time scales are of the order 
of nanoseconds, 10-100 psec, and 100 fsec-1 psec. If we consider a proton at the same location as the chromophore 
and described by a continuous co-ordinate x, then the time-dependent interaction energy with the reaction field is 
ex(t)R(t), where e is the proton charge. For the chromophore, modelled as a two-level system the interaction energy 
is AfiR(t)(P e — P g ){t) where the last factor represents the relative occupation probability of the excited and ground 
states of the chromophore. Then the spectral density, J(u>) = MRe 7(w)w, relevant to the Caldeira-Leggett model 
can be related to the spectral density J c {oj) associated with a spin boson model for the chromophore, @ 

JH = (i) Jc{uj) - (M) 

Comparison with ultra-fast laser spectroscopy experiments, with molecular dynamics simulations, and continuum di- 
electric models shows that@ the high-frequency (w > 10 cm -1 ) part of the chromophore spectral density is dominated 
by the bulk water. 

It should be noted that the femtosecond laser spectroscopy experiments used to extract spectral densities do have 
limited time resolution (~ 10-100 fsec) and so only give information about the spectral density for u> < 500 cm -1 . 
Furthermore, quantum Kra mers theory requires a knowledge of 7(w) and "f(z) at frequencies of order loq and uib- 

Fig. 1 of Lang et al. [l09j | shows the frequency dependence of a chromophore- w ater spectral density for frequencies 
up to about ui ~ 3000 cm -1 , calculated by the method of Song and Chandler. [llOj It has peaks of magnitude of 
the order of a few hundred cm -1 , at frequencies of about 180, 600, 1800, and 3200 cm -1 . There is a substantial 
contribution from the librational band of water centred at 600 cm . 

If we combine (|2Tj) with (|54|) we find that the curvature K e can be related to the reorganisation energy Er associated 
with the excited state of the chromophore 0, 



and 



2 



*<*^dbGy"** (56) 

Given typical values of Er of ~ 1000 cm -1 6] and A/i ~ 5 D, for chromophores in proteins ((he/ A/i) 2 /M ~ 30 cm -1 
for protons and A/i = 1 D) (|56|) then gives an estimated upper bound 
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Hence, we see that for typical protein and solvent environments we expect to be in the weak friction limit for proton 
transfer reactions at room temperature, since according to the discussion around Eqn. (|46[) only values of z ~ u>b and 
larger are relevant. 

2. Infra-red spectroscopy 

Infra-red (IR) spectroscopy provides a means to measure the frequencies, damping, and decoherence of vibrational 
modes in a molecule. Rece nt ad vances in femtosecond two-dimen siona l IR spectroscopy has yielded such inform ation 
for several specific systems. [lll| including O-H stretches in water, pjjj N-H stretches in Watson-Crick base pairs |l!3j . 
It is generally found that hydrogen bonding leads to broad IR spectra. [102L Ill2l Ill4| 

3. Molecular dynamics simulations 

Stimulated by recent time-resolved infra-red spectroscopy measurements of vibrational energy relaxation and de- 
phasing a number of molecular dynamics studies have been made for specific vibrations in biomolcculcs. The difficulties 
associated with extracti ng th e vibrational energy relaxation rate, 1/T± from molecular dynamics simulations has re- 
cently been summarised. [115J If the vibration frequency is in the classical regime (ujq < ksT ~ 200 cm" 1 ) then the 
Landau- Teller equation allows one to extract T\ from the classical forc e-forc e correlation function. However, most 
modes of interest are not in this classical regime. Fujisaki and Straub jll5l | considered the specific case of a C-D 
vibration in cytochrome-c surrounded by water. The experimental value of T\ for this mode is about one psec. They 
found that the value of T\ found in the simulations could vary by as much as two orders of magnitude with only a ten 
per cent change in the bond force constant. Time resolved infra-red spectroscopy experiments show that the C=0 
stretch of the peptid e bo nd for a wide range of proteins has a relaxation time of about one psec. A recent molecular 
dynamics simulation [ll6j | yielded values that were two orders of magnitude larg er than this. 

A recent combined molecular dynamics-quantum mechanics calculation [1171 ] calculated the IR spectrum for the 
water networks in (the proton pump) bacteriorhodopsin. They found broad continuum bands, around 1800 and 2700 
cnr 1 , and associated them with the solvated Zundel complex (H5OJ) and Eigen complex (HaO + ). A time-resolved 
fourier-transform infrared spectroscopy experimen t fou nd that the precise arrangement of the water molecules within 
bacteriorhodopsin was crucial to proton transfer. [l!8j | Given that many proton transfer reactions in enzymes also 
involve water molecules inside the protein, such hydrogen bonding networks may also be a significant source of 
friction for proton transfer. 

Moritsugu and Smith performed molecular dynamics simulations of myoglobin bot h in w ater and in vacuum. They 
used a Langevin model to describe the dynamics of the different vibrational models, jl lOj ] The frictional damping of 
different vibrational modes of the protein was found to be proportional to the accessible surface area of the mode, 
confirming the importance of the solvent which we have stressed here. At 300 K, they found that for modes with u> 
in the range 100-400 cm -1 , the friction could be fit to 

Re j(u) = A7 + Aw (58) 

with A7 ~ 20 cm -1 , and A ~ 0.38 and was temperature independent between 120 and 300 K. 

4- Dielectric continuum models 

Continuum models 0, IllOL Il2dl | allow one to express the spectral density in terms of the frequency dependent 
(complex) dielectric function e s (u>) of the solvent. If the proton is at the centre of a spherical cavity of radius a inside 
the water, and undergoes displacements much less than the radius, then using (|54|) and results for chromophores then 

Re 7 M = - e \ Im ( ^\~ ec ) (59) 
zneoa^Mui \2e s (uj) + e c J 

and e c is the (static) dielectric constant of the cavity, which can be approximated as the local dielectric constant of 
the protein environment surrounding the proton. 

Measurements of the frequency dependent dielectric constant of water e s (w) in the r ange 1-200 cm , have been fit 
to a form involving three Debye relaxation terms and one damped resonant term [l2l| 



(60) 
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where r, is the Debye relaxation time of the relevant component. For water at 298 K, the co-efficients are Ae^ 
(i = 1,2,3,4) = 71.5, 2.8, 1.6, and 0.92, respectively. The correspondin g rela xation times are Tj = 8.3, 1.0, 0.1, and 
0.025 psec. The resonant frequency is = 175 cm -1 . Roy and Bagchi |l22j gave a resonant frequency of u>i = 200 
cm -1 and a damping constant such that W4T4 = 2. They also calculated the frequency dependent friction for outer 
sphere electron tr ansfe r reactions in water out to lot — 5 where r = 0.1 psec. 

Figure 5 of Ref. Il2ll shows that for u> ~ 100 cm -1 that Imefw) ~ 2 (see also Fig. 18 in Ref. Il23l ) which is an order 
of magnitude larger than the contribution from the slowest relaxation (i = 1). 



A parametrisation of the higher frequency part of the dielectric function has been given. 124j There are features at 
frequencies of about 180, 600, 1800, and 3200 cm -1 . The first can be assigned to hindered translation of the hydrogen- 
bonded network (the O- • • O stretch of the O-H - • • O of the water hydrogen bonds). Hindered rotation (libration) is 
the origin of the second feature. Fig. 3 in Ref. Il20l shows a plot of the frequency dependence of the right hand site 
of ([55]) . Again, there is a substantial contribution from the librational band in the range 600-800 cm -1 . 

We can model a peak at a frequency Lo r in the spectral density by ([23]) with the corresponding friction kernel, (|24[) . 
If z, to r 3> r then 

7(z)~^. (61) 

and for the typical values of j r ~ T ~ 100's cm -1 , ujb ~ v — 2nkBT ~ 1000 cm -1 this gives *f(v) <C u which justifies 
using the weak-friction limit in reaction rate theory. 



VII. QUANTUM TRANSITION STATE THEORY DESCRIBES THE EXPERIMENTAL DATA 

Quantum Transition State Theory predicts that for T > To the H/D KIE is given by 

fcg _ /g smhjfeVgfcgT) sm(hw b /2V2k B T) 
k D sinli(fiw /2 v / 2fc B T) sm(haj b /2k B T) 

Note that this expression depends only on two parameters, loq and u>b- Figures [3] and |4] show the fits of Eq. (|62|) 
to experimental data for different enzymes that show kinetic anomalies typical of systems that have been argued to 
exhibit tunneling. The first thing to notice is that the quantum transition state theory result can reproduce these 
experimental results. Secondly, the values of loq and loj, obtained from the fits are comparable to typical C-H stretch 
frequencies and to barrier frequencies that have been obtained from quantum chemistry calculations (Table [TTTJ) . The 
values of LOb obtained for both fits imply a crossover temperature To below room temperature, indicating that our 
description is self-consistent in that we are in the temperature regime above the crossover. 



VIII. CONCLUSIONS 



The path integral approach we have used has several benefits for elucidating questions concerning the role of 
quantum tunneling in hydrogen transfer reactions in enzymes. It allows a full quantum mechanical treatment of the 
role of both temperature and the environment of the active site. Quantum tunneling is described by the instantons 
(bounce solution), which are periodic solutions to the semi-classical equations of motion in imaginary time. An 
important result is that these solutions only exist below some temperature, To which is determined by the curvature 
of the top of the barrier. Above this temperature the only role of quantum effects concerns quantum fluctuations about 
the transition state. We have shown that for two enzymes a quantitative description of the temperature dependence 
of kinetic isotope effects is possible in terms of such a quantum transition state theory. This suggests that quantum 
tunneling does not play as significant role in hydrogen transfer as is often claimed. 
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FIG. 3: Comparison of the measured temperature dependence of the kinetic isotope effect [43l] for the enzyme methylmalonyl- 
CoA mutase and several synthetic molecules with similar H-atom abstraction reactions (NpCb, AdoCbl and 8-MeOAdoCbl) 
with Quantum Transition State Theory at temperatures above which tunneling is possible. These systems all have KIE's and 
Arrhenius parameters that differ by factors of 5-10 from the semi-classical values traditionally claimed to hold in the absence 
of tunneling (see Table I). The use of the synthetic molecules allows coverage of a much wider temperature range (from 10 to 
120° C) than possible with enzymes. The solid line is a fit of Eqn. (|62|) to the experimental data with two free parameters, 
the barrier frequency, u)b, and the oscillation frequency in the reactant well, u>o- The value obtained for u>q is comparable with 
typical C-H stretch frequencies. The value obtained for uib is comparable with estimates from quantum chemistry calculations 
(see Table ITTTjl of similar reactions. The value of the crossover temperature To ~ 250 K = — 20°C implied by the fitted value 
of u>b and Eqn. |[8]| is consistent with the domain of validity of (|62|l (i.e., T > To). This Figure shows that it is not necessary 
to invoke quantum tunneling to obtain a quantitative description of the experimental data for this enzyme and its synthetic 
analogues. 



APPENDIX A: THE SWAIN-SCHAAD EXPONENT IS AN UNRELIABLE CRITERIA FOR QUANTUM 

TUNNELING 



A parameter that has been used to quantify the effect of isotope substitution is the Swain-Schaad exponent, which 
is defined as 

_ ln(WM , , 

a -Hk D /k T y (A1) 

This parameter is used because in the standar d a pproach to reaction rate theory the Swain-Schaad exponent is 
independent of temperature and system specifics [1251 ] . It can be written as 

_ \n(A H /A T ) - [E P - E T ] /k B T 
a \n{A D /A T )-[E D -E T ]/k B r [ ' 

In the limit of low temperatures and where Ah /At and A o j At — > 1 this simplifies to 

\n(k H /k T ) E P -E T 



\n{k D /k T ) E D - E T 



3.3 (A3) 



From the standard semi-classical rate theory the ratio of the differences in activation energies of the different 
isotopes is a constant that does not depend on the temperature or system specifics. In this case the value that a 
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FIG. 4: Comparison of the temperature dependence of the kinetic isotope effect (for tritium substitution) measured for fiavoen- 
zyme monoamine oxidase [2f| with Quantum Transition State Theory at temperatures above which tunneling is possible. The 
solid line is a fit of Eqn. (|47[) to the experimental data with two free parameters, the barrier frequency, ut, and the oscillation 
frequency in the reactant well, too- The value obtained for ujb is comparable with estimates from quantum chemistry calculations 
(see Table HIT)) for an amineoxidase enzyme which give a barrier frequency around 1000 cm -1 [9q |. The value of 2100 cm -1 for 
the reactant well oscillation frequency is about two-thirds of typical carbon hydrogen stretch frequencies, but is comparable to 
values of the stretch frequency calculated near transition states [55l. l59j|. The value of the crossover temperature To ~ 240 K 
implied by the fitted value of uib and Eqn. ([8J is consistent with the domain of validity of (I62|l . This Figure shows that it is 
not necessary to invoke quantum tunneling to obtain a quantitative description of the experimental data for this enzyme. 

takes will depend on the prefactor values and also on temperature. This subtlety has sometimes been overlooked 
in the experimental literature. There are a number of cases where a discrepancy between the theoretical value of 
(Ep — Et)/{Ed — Et) and the experimentally determined value oi\n(kH / kx) / ^{ko / kr) have formed the basis for the 
hypothesis that tunneling is occurring, even though all the terms in the above expression are comparable [U HH, EH- 
It is for this reason that in this paper we do not use the Swain-Schaad exponent as a criteria for the presence of 
tunneling. 



APPENDIX B: EFFECT OF THE ENVIRONMENT ON THE CROSSOVER TEMPERATURE 

In the case of a Lorentzian spectral density equation for /i takes the form 

M - u b + ■ — = 0. (Bl) 

To varies as a function of the friction strength for different bath response frequencies, u)d- Fig. [5] shows how the 
positive root of Eq. (|B1[) changes as a function of the scaled friction and bath response frequency. For all values 
of lvd an increase in the strength of the damping, 7, reduces the effective barrier frequency and hence the crossover 
temperature. 
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FIG. 5: The crossover temperature for a parabolic barrier as a function of the friction strength for a Lorentzian spectral 
density. The friction strength is the dimensionless parameter 7/0^6. The different plots show how the relative size of the 
response frequency of the bath, u)d, and the barrier frequency, u>t change the influence that friction has on the crossover 
temperature. These show that when lub 3> <^>b the crossover temperature becomes most sensitive to the friction strength. Then 
the bath is able to respond on a timescale which is much faster than all of the other relevant timescales. When the response 
frequency of the bath is much less than the barrier frequency the crossover temperature is only weakly suppressed as the friction 
strength increases. 



G. R. Fleming and G. D. Scholes, Nature 431, 256 (2004). 

D. Abbott, P. Davies, and A. Pati, eds., Quantum Aspects of Life (Imperial College Press, London, 2008). 
V. Helms, Curr. Opin. Struct. Biol. 12, 169 (2002). 

I. Tejero, N. Gonzalez-Garcia, A. Gonzalez-Lafont, and J. M. Lluch, J. Am. Chem. Soc. 129, 5846 (2007). 

J. Gilmore and R. McKenzie, J. Phys. Chem. A 112, 2162 (2008). 

A. Nagy, V. Prokhorenko, and R. Miller, Curr. Opin. Struct. Bio. 16, 654 (2006). 

J. Reimers, L. McKemmish, A. Mark, R. McKenzie, and N. Hush, Proc. Nat. Acad. Sci. (USA) 106, 4219 (2009). 
L. McKemmish, J. Reimers, A. Mark, R. McKenzie, and N. Hush, Phys. Rev. E 80, 021912 (2009). 

G. Engcl, T. R. Calhoun, E. L. Read, T.-K. Ahn, T. Mancal, Y.-C. Cheng, R. E. Blankenship, and G. R. Fleming, Nature 
446, 782 (2007). 

R. Allemann and N. Scrutton, eds., Quantum Tunneling in Enzyme-Catalyzed Reactions (Royal Society of Chemistry, 
London, 2009). 

U. Weiss, Quantum Dissipative Systems (World Scientific, Singapore, 1999), 2nd ed. 
M. Schlosshauer, Rev. Mod. Phys. 76, 1267 (2004). 

C. C. Moser, C. C. Page, and P. L. Dutton, Phil. Trans. R. Soc. B 361, 1295 (2006). 

D. Devault, Photosyn. Res. 22, 5 (1989). 

N. Alberding, R. H. Austin, K. W. Beeson, S. S. Chan, L. Eisenstein, H. Frauenfelder, and T. M. Nordlund, Science 192, 
1002 (1976). 

C. Ober, M. Burkardt, H. Winkler, A. X. Trautwein, A. A. Zharikov, S. F. Fischer, and F. Parak, Eur. Biophys. J. 26, 
227 (1997). 

R. Marcus, J. Chem. Phys. 125, 194504 (2006). 

R. A. Marcus, Phil. Trans. R. Soc. B 361, 1445 (2006). 

P. L. Dutton, A. Munro, N. S. Scrutton, and M. J. Sutcliffe, Phil. Trans. R. Soc. B 361, 1293 (2006). 

J. P. Klinman, Phil. Trans. R. Soc. B 361, 1323 (2006). 

Y. Cha, C. J. Murray, and J. P. Klinman, Science 243, 1325 (1989). 



20 



[23] K. L. Grant and J. P. Klinman, Biochemistry 28, 6597 (1989). 

[24] B. J. Bahnson, D. Park, K. Kim, B. V. Plapp, and J. P. Klinman, Biochemistry 32, 5503 (1993). 

[25] T. Jonsson, D. E. Edmondson, and J. P. Klinman, Biochemistry 33, 14871 (1994). 

[26] T. Jonsson, M. H. Glickman, S. J. Sun, and J. P. Klinman, J. Am. Chem. Soc. 118, 10319 (1996). 

[27] J. C. Nesheim and J. D. Lipscomb, Biochemistry 35, 10240 (1996). 

[28] M. M. Whittaker, D. P. Ballou, and J. W. Whittaker, Biochemistry 37, 8426 (1998). 

[29] J. Basran, M. J. Sutcliffe, and N. S. Scrutton, Biochemistry 38, 3218 (1999). 

[30] A. Kohen, R. Cannio, S. Bartolucci, and J. P. Klinman, Nature 399, 496 (1999). 

[31] S. Chowdhury and R. Banerjee, J. Am. Chem. Soc. 122, 5417 (2000). 

[32] R. Harris, R. Meskys, M. J. Sutcliffe, and N. S. Scrutton, Biochemistry 39, 1189 (2000). 

[33] J. L. Abad, F. Camps, and G. Fabrias, Angew. Chem. Int. Ed. 39, 3279 (2000). 

[34] S. L. Seymour and J. P. Klinman, Biochemistry 41, 8747 (2002). 

[35] N. Angrawal, B. Y. Hong, C. Mihai, and A. Kohen, Biochemistry 43, 1998 (2004). 

[36] K. Remmers, W. L. Meerts, and I. Ozier, J. Chem. Phys. 112, 10890 (2000). 

[37] P. R. Schreiner, H. P. Reisenauer, F. C. P. IV, A. C. Simmonett, W. D. Allen, E. Mtyus, and A. G. Csszr, Nature 453, 
906 (2008). 

[38] V. A. Benderskii, V. I. Goldanskii, and D. E. Makarov, Phys. Rep. 233, 195 (1993). 

[39] V. A. Benderskii, D. E. Makarov, and C. A. Wight, Chemical Dynamics at Low Temperatures, vol. 88 of Advances in 

Chemical Physics (Wiley, New York, 1994) . 

[40] M. Garcia- Viloca, J. Gao, M. Karplus, and D. G. Truhlar, Science 303, 186 (2004). 

[41] J. Villa and A. Warshel, J. Phys. Chem. B 105, 7887 (2004). 

[42] A. Kohen and J. P. Klinman, Acc. Chem. Res. 31, 397 (1998). 

[43] K. M. Doll and R. Finke, Inorg. Chem. 42, 4849 (2003). 

[44] K. M. Doll, B. R. Bender, and R. G. Finke, J. Am. Chem. Soc. 125, 10877 (2003). 

[45] S. Caratzoulas, J. S. Mincer, and S. D. Schwartz, J. Am. Chem. Soc. 124, 3270 (2002). 

[46] P. Bertrand, J. Biol. Inorg. Chem. 9, 2 (2004). 

[47] D. Antoniou, J. Basner, S. Nunez, and S. Schwartz, Chem. Rev. 106, 3170 (2006). 

[48] J. Pu, J. Gao, and D. Truhlar, Chem. Rev. 106, 3140 (2006). 

[49] K. Peters, Acc. Chem. Res. 42, 89 (2009). 

[50] M. J. Knapp, K. Rickert, and J. P. Klinman, J. Am. Chem. Soc. 124, 3865 (2002). 

[51] M. P. Meyer and J. P. Klinman, Chem. Phys. 319, 283 (2005). 

[52] E. Hatcher, A. Soudackov, and S. Hammes-Schiffer, J. Phys. Chem. B 109, 18565 (2005). 

[53] E. Hatcher, A. V. Soudackov, and S. Hammes-Schiffer, J. Am. Chem. Soc. 129, 187 (2007). 

[54] W. Siebrand and Z. Smedarchina, J. Phys. Chem. B 108, 4185 (2004). 

[55] Q. Cui and M. Karplus, J. Am. Chem. Soc. 124, 3093 (2002). 

[56] M. Wang, Z. Lu, and W. Yang, J. Chem. Phys. 124, 124516 (2006). 

[57] Q. Cui, M. Elstner, and M. Karplus, J. Phys. Chem. B 106, 2721 (2002). 

[58] C. Alhambra, J. C. Corchado, M. L. Sanchez, J. Gao, and D. G. Truhlar, J. Am. Chem. Soc. 122, 8197 (2000). 

[59] C. Alhambra, M. L. Sanchez, J. C. Corchado, J. Gao, and D. G. Truhlar, Chem. Phys. Lett. 355, 388 (2002). 

[60] P. F. Faulder, G. Tresadern, K. K. Chohan, N. S. Scrutton, M. J. Sutcliffe, I. H. Hillier, and N. A. Burton, J. Am. Chem. 
Soc. 123, 8604 (2001). 

[61] G. Monard, X. Prat-Resina, A. Gonzalez-Lafont, and J. M. Lluch, Int. J. Quantum Chem. 93, 229 (2003). 

[62] S. Hammes-Schiffer, Curr. Opin. Struct. Biol. 14, 192 (2004). 

[63] G. Tresadern, J. P. McNamara, M. Mohr, H. Wang, N. A. Burton, and I. H. Hillier, Chem. Phys. Lett. 358, 489 (2002). 

[64] L. Masgrau, A. Roujeinikova, L. O. Johannissen, P. Hothi, J. Basran, K. E. Ranaghan, A. J. Mulholland, M. J. Sutcliffe, 

N. S. Scrutton, and D. Leys, Science 312, 237 (2006). 

[65] M. H. M. Olsson, P. E. M. Siegbahn, and A. Warshel, J. Am. Chem. Soc. 126, 2820 (2004). 

[66] D. Antoniou, S. Caratzoulas, C. Kalyanaraman, J. S. Mincer, and S. D. Schwartz, Eur. J. Biochem. 269, 3103 (2002). 

[67] A. Caldeira and A. Leggett, Ann. Phys. 149, 374 (1983). 

[68] P. Hanggi, P. Talkner, and M. Borkovec, Rev. Mod. Phys. 62, 251 (1990). 

[69] H. Eyring, J. Chem. Phys. 3, 107 (1935). 

[70] Y. H. Kim and M. M. Kreevoy, J. Am. Chem. Soc. 114, 7116 (1992). 

[71] R. Bell, The Tunnel Effect in Chemistry (Chapman and Hall, London, 1980). 

[72] A. Kohen and H. H. Limbach, eds., Isotope Effects in Chemistry and Biology (Taylor & Francis, Boca Raton, 2006), 

chap. 28, pp. 743-764. 

[73] E. Wigner, Z. Phys. Chem. Abt. B. 19, 203 (1932). 

[74] M. P. Meyer, D. R. Tomchick, and J. P. Klinman, Proc. Nat. Acad. Sci. (USA) 105, 1146 (2008). 

[75] J. Basran, S. Patel, M. J. Sutcliffe, and N. S. Scrutton, J. Biol. Chem. 276, 6234 (2001). 

[76] J. Basran, M. J. Sutcliffe, and N. S. Scrutton, J. Biol. Chem. 276, 24581 (2001). 

[77] W. A. Francisco, M. J. Knapp, N. J. Blackburn, and J. P. Klinman, J. Am. Chem. Soc. 124, 8194 (2002). 

[78] R. Sikorski, L. Wang, K. A. Markham, P. T. R. Rajagopalan, S. J. Benkovic, and A. Kohen, J. Am. Chem. Soc. 126, 
4778 (2004). 

[79] G. Maglia and R. Allemann, J. Am. Chem. Soc. 125, 13372 (2003). 

[80] N. Fujisaki, A. Ruf, and T. Gaumann, J. Am. Chem. Soc. 107, 1605 (1985). 



21 



N. Fujisaki, A. Ruf, and T. Gaumann, J. Chem. Phys. 80, 2570 (1983). 

J. Braun, R. Schwesinger, P. G. Williams, H. Morimoto, D. E. Wemmer, and H. H. Limbach, J. Am. Chem. Soc. 118, 
11101 (1996). 

E. F. Caldin and S. Mateo, J. Chem. Soc., Chem. Commun. 22, 854 (1973). 
R. Zwanzig, J. Stat. Phys. 9, 215 (1973). 

J. Ankerhold, Quantum tunneling in complex systems: the semi-classical approach (Springer, Berlin, 2007). 
H. Grabert, P. Olschowski, and U. Weiss, Phys. Rev. B 36, 1931 (1987). 
W. H. Miller, J. Chem. Phys. 62, 1899 (1975). 

S. Coleman, Aspects of Symmetry (Cambridge University Press, 1995). 
W. H. Miller, Adv. Chem. Phys. 30, 77 (1975). 

E. Freidkin, P. Riseborough, and P. Hanggi, Z. Phys. B 64, 237 (1986). 

P. S. Riseborough, P. Hanggi, and E. Freidkin, Phys. Rev. A 32, 489 (1985). 

P. G. Wolynes, Phys. Rev. Lett. 47, 968 (1981). 

P. Hanggi, Z. Phys. B: Condens. Matter 68, 181 (1987). 

F. Settle, ed., Handbook of Instrumental Techniques for Analytical Chemistry (Prentice-Hall, Upper Saddle River, 1997), 
chap. 15, p. 266. 

C. Rotting and K. Gerwert, Chem. Phys. 307, 227 (2004). 

S. S. Erdcm, O. Karahan, I. Yildiz, and K. Yelekci, Org. Biomol. Chem. 4, 646 (2006). 

S. Nunez, G. Tresadern, I. Hillier, and N. Burton, Phil. Trans. R. Soc. B 361, 1387 (2006). 

L. Masgrau, K. E. Ranaghan, N. S. Scrutton, A. J. Mulholland, and M. J. Sutcliffe, J. Phys. Chem. B 111, 3032 (2007). 

S. J. Titmuss, P. L. Cummins, A. A. Bliznyuk, A. P. Rendell, and J. E. Gready, Chem. Phys. Lett. 320, 169 (2000). 

A. Dybala-Defratyka, P. Paneth, R. Banerjee, and D. G. Truhlar, Proc. Natl. Acad. Sci. 104, 10774 (2007). 

V. May and O. Kuhn, Charge and Energy Transfer Dynamics in Molecular Systems (Wiley- VCH, Berlin, 2000). 

K. Giese, M. Petkovic, H. Naundorf, and O. Kuhn, Phys. Rep. 430, 79 (2006). 

N. Rosch and M. Ratner, J. Chem. Phys. 61, 3344 (1974). 

J. Eaves, A. Tokmakoff, and P. Geissler, J. Phys. Chem. A 109, 9424 (2005). 

C. Fecko, J. Eaves, J. Loparo, A. Tokmakoff, and P. Geissler, Science 301, 1698 (2003). 

C. Fecko, J. Loparo, S. Roberts, and A. Tokmakoff, J. Chem. Phys. 122, 054506 (2005). 

O. Kuhn and H. Naundorf, Phys. Chem. Chem. Phys. 5, 79 (2003). 

A. Nitzan, Chemical dynamics in condensed phases (Oxford University Press, Oxford, 2006). 
M. J. Lang, X. J. Jordanides, X. Song, and G. R. Fleming, J. Chem. Phys. 110, 5884 (1999). 
X. Song and D. Chandler, J. Chem. Phys. 108, 2594 (1998). 
M. Cho, Chem. Rev. 108, 1331 (2008). 

S. Woutersen and H. Bakker, Phys. Rev. Lett. 83, 2077 (1999). 
S. Woutersen and G. Cristalli, J. Chem. Phys. 121, 5381 (2004). 
T. Yamashita and K. Takatsuka, J. Chem. Phys. 126, 074304 (2007). 
H. Fujisaki and J. Straub, Proc. Nat. Acad. Sci. (USA) 102 (2005). 
P. Nguyen and G. Stock, J. Chem. Phys. 119, 11350 (2003). 

G. Mathias and D. Marx, Proc. Nat. Acad. Sci. (USA) 104 (2007). 
F. Garczarek and K. Gerwert, Nature 439, 109 (2006). 

K. Moritsugu and J. Smith, J. Phys. Chem. B 109, 12182 (2005). 

C.-P. Hsu, X. Song, and R. A. Marcus, J. Phys. Chem. B 101, 2546 (1997). 

K. Vij, D. Simpson, and O. Panarina, J. Mol. Liq. 112, 125 (2004). 

S. Roy and B. Bagchi, J. Chem. Phys. 99, 9938 (1993). 

H. Hasted, Water: a comprehensive treatise (Plenum, New York, 1972), chap. 7. 
N. Nandi, S. Roy, and B. Bagchi, J. Chem. Phys. 102, 1390 (1995). 

C. G. Swain, E. C. Stivers, J. F. Reuwer, and L. J. Schaad, J. Am. Chem. Soc. 80, 5885 (1958). 



